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Abstract 

■^h ■ We present an accurate ab initio method of calculating isotope shifts and relativistic shifts in 

£>v atomic spectra. We test the method on neutral carbon and three carbon ions. The relativistic 



shift of carbon lines may allow them to be included in analyses of quasar absorption spectra 
that seek to measure possible variations in the fine structure constant a over the lifetime of the 
Universe. Carbon isotope shifts can be used to measure isotope abundances in gas clouds: isotope 
abundances are potentially an important source of systematic error in the a-variation studies. 
These abundances are also needed to study nuclear reactions in stars and supernovae, and test 
! models of chemical evolution of the Universe. 
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I. INTRODUCTION 



Recent studies of quasar absorption systems have suggested that the fine structure con- 
stant, a = e 2 /hc, may have been smaller in the early universe than it is today on Earth 
IH • Other groups using different telescopes have shown zero variation (H El, 0| . All 
of these studies use the "many multiplet" method [8|, which relies on the comparison of the 
wavelengths of many transitions in many ions to enhance the size of the effects and remove 
sources of systematic error. While this method offers an order-of-magnitude improvement 
in sensitivity over the previously used alkali-doublet method, a potential systematic error is 
introduced related to the isotope abundances of the absorbers. If the isotope abundances of 
the absorbers differ from the terrestrial abundances then there can be spurious shifts in the 
measured wavelengths which may be incorrectly interpreted as variation of a. 

This problem can be resolved if both the relativistic shift and isotope shift of each tran- 
sition is known, by using particular combinations of the transition frequencies as probes 
which are insensitive to either a- variation or isotopic abundances 0. Changes in the iso- 
tope abundances and the fine-structure constant can then be measured directly. 

The measured isotopic abundance of carbon can be used to test models of chemical 
evolution of the Universe. Of particular importance are some chemical evolution models 
with enhanced populations of intermediate-mass stars that serve as factories for heavy Mg 
isotopes. The changes in the relative abundances of Mg isotopes could account for much of 



the evidence for variation in a at relatively low redshifts (z < 2) |1QL lll|. However, these 
models also overproduce nitrogen, violating observed abundance trends in high-z damped 
Lyman-a systems. Furthermore, it has been shown that such models would also increase 
the ratio of 13 C to 12 C [12] . With the isotope shifts calculated in this paper it is possible to 
measure this ratio, and hence provide a diagnostic of these non-standard chemical evolution 
models. 

We also have one more reason to study carbon: it is a well studied atom, and we can 
compare the results of our method with those of other theoretical analyses, as well as a 
few experiments (unfortunately, isotope shifts for the majority of lines used in astrophysical 
applications have not been measured). In particular, much progress has been made to cal- 
culate isotope shifts using the multiconfigurational Hartree-Fock (MCHF) and configuration 
interaction (CI) approach [IH, EH EH EH • I n Sections HP and ITTT1 we present our method of 
calculating the isotope shift. It uses the finite-field method to reduce the problem of isotope 
shift (or relativistic shift) to that of an energy calculation. The transition energies are cal- 
culated using a combination of the configuration interaction and many-body perturbation 
theory (MBPT) approaches, following the works of Refs. E3 an d E3- 

In Sec. II VI we show that our CI + MBPT method has accuracy comparable to that of the 
MCHF calculations, and are accurate enough to be used to measure isotope abundances. 
We believe that our method may have some advantages over the MCHF calculations; in 
particular, it is more readily applicable to heavier atoms (where we need to calculate isotope 
shifts for the study of the a variation and isotopic evolution) and has already been shown to 
be more accurate in the case of Mg I ^| . We also use this method to calculate dependence 



of the carbon transition frequencies on a, which is needed to study a variation. 
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II. METHOD 



Using many-body perturbation theory in the residual Coulomb operator and specific 
mass shift (SMS) operator to calculate isotope shift shows poor convergence. Therefore, 
we are looking for an "all order" method of calculation. The finite-field scaling method is 
used, reducing the task to an energy calculation, and including the SMS in all parts of the 
calculation. A similar idea is used to calculate relativistic shift. 



A. Relativistic shift 

To measure a in the distant past, we compare the frequency of transitions in the labo- 
ratory, uq, with those measured in the quasar absorption spectra, lu. The difference can be 
expressed as 

to = lu + qx, (1) 

where x = (a/ do) 2 — 1 and cto is the laboratory value of the fine structure constant. We 
vary a directly in computer codes, and extract the relativistic shift q from the calculated 
spectrum as 

q = ~r • 2 

dx x =o 

Thus we have reduced the problem to an accurate numerical calculation of the energy, and 
hence u. 



B. Isotope shift 

The isotope shifts of atomic transition frequencies come from two sources: the finite size 
of the nuclear charge distribution (the "volume" or "field" shift), and the finite mass of 
the nucleus (see, e.g., Ref. [ZD]). In atoms as light as carbon, the field shift is negligible in 
comparison to the mass shift: we will not consider it further. 

Because a real nucleus has a finite mass, there is a recoil effect from the movements of 
the electrons. The energy shift due to the recoil of the nucleus of mass M is given by 

2M 2M Yi 7 

This "mass shift" is traditionally divided into the normal mass shift (NMS) and the specific 
mass shift (SMS), given by the first and second terms of Eq. (j3J), respectively. The normal 
mass shift is easily calculated from the transition frequency; the specific mass shift is difficult 
to evaluate accurately. 

The shift in energy of any transition in an isotope with mass number A' with respect to 
an isotope with mass number A can be expressed as 

bv A!,A = V A< _ y A = (a , nms + hus) (1 - I) , (4) 
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where the normal mass shift constant is 



fcNMS = -^3. (5) 

The value 1823 refers to the ratio of the atomic mass unit to the electron mass. 

To calculate k SMS we include a scaled specific-mass- shift operator directly into our energy 
calculation from the very beginning. We add the two-body SMS operator to the Coulomb 
potential Q = 1/ \ri — r 2 | + Api • p 2 . The operator Q replaces the Coulomb operator 
everywhere that it appears in an energy calculation. We recover the specific mass shift 
constant as 

&SMS = -77 (6) 

CLA A=0 

The operator Q has the same symmetry and structure as the Coulomb operator (see Ap- 
pendix EJ. 



III. ENERGY CALCULATION 

To calculate the energies we first solve the Dirac-Fock equations for the core and valence 
electrons, and we generate a basis set that includes the core and valence orbitals and a 
number of virtual orbitals. We then calculate the energy levels using a combination of CI 
and MBPT for many- valence-electron atoms as was done for Mgl in Ref. 0]. In this section 
we outline the procedure for the combined CI and MBPT method; it generally follows the 



work of Ref. [17(. Note that for single- valence-electron atoms, where the CI procedure is 



unnecessary, the method reduces to the addition of core-correlation effects to the Dirac-Fock 



energy using MBPT, as was shown to be highly successful for calculating SMS in Ref. |21 
Atomic units (h = m e = e = 1) are used throughout this paper. 



A. Single particle basis 

We firstly solve the Dirac-Fock equations for one-particle wavefunctions of the open-shell 
core \m) 

h DF \m) = e m \m) , (7) 

where 

h DF = ccx-p+(p-l)m e c 2 -- + V N ^(r) (8) 

r 

where V Ndf is the potential (both direct and exchange) of the Adf electrons included in 
the self-consistent Hartree-Fock procedure. Note that this is not necessarily the number of 
electrons in the closed-shell core, in fact N com < Adf < A, where N is the total number of 
electrons. For the purposes of the CI calculation there are N — N com valence electrons. 

We need to generate a basis set, |z) that includes the core and valence states and a 
large number of virtual states. In this paper we have used a B-spline basis set, formed by 
diagonalizing the Dirac-Fock operator on the basis set of B-splines and excluding orbitals 
with high energy. This basis has been shown to be effective in energy calculations using this 
method of CI and MBPT Q. 
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B. Configuration interaction method 

The many-electron Hilbert space is separated into subspaces P and Q. P corresponds to 
the frozen-core approximation; Q is complimentary to it and includes all states with core 
excitations. Using Slater determinants (denoted with capital letters) \I) of the single-particle 
functions \i) as a basis set in the many-electron space, we can define projection operators 
for P and Q by 

7>=EI J > (!\ (9) 

leP 

V + Q = l (10) 

Determinants that have all core states fully occupied are in the P subspace; all others are 
in the subspace Q. 

In the CI method, the many-electron wavefunction is expressed as a linear combination 
of Slater determinants from the subspace P: 

^ = T,Ci\i), (ii) 
i&p 

where the Cj are obtained from the matrix eigenvalue problem 

E HuCj = ECi. (12) 
JeP 

In the frozen-core approximation, the determinants \I) need include only the valence elec- 
trons. Although P is infinite-dimensional, we can use a finite-dimensional model subspace 
by specifying the set of allowed configurations for the valence electrons, for example by re- 
stricting the set of single particle states. The restrictions we use are different for each ion, 
and are expressed more fully in the relevant parts of Section IIVI We will not distinguish 
here between the finite and infinite subspaces. 

The Hamiltonian for the CI problem is a projection of the exact Hamiltonian 7i onto the 
model subspace. The core is frozen in the P subspace, so our projected Hamiltonian is 

vHv = E cove + e h ? l + E — ( 13 ) 

i>N COIC j>i>N COIC V 

where E core is the total energy of the N core core electrons and the single particle operator 

h CI = c a ■ p + (B - l)m e c 2 - - + ^ Ncoro (14) 

r 

acts only on the valence electrons. Using the operator ()13|) in Eq. (|T2*|) corresponds to the 
pure CI method in the frozen-core approximation. 



C. Exact Hamiltonian expansion 

We wish to write the exact equivalent of 7i in the subspace P. The "Feshbach operator" 
yields the exact energy when operating on the model function \l/p = '. Following Lindgren 
and Morrison j^, we start from the many-body Schrodinger equation in the form 

H(V+Q)V = EV (15) 
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and operate from the left with V and Q, respectively to obtain 



VHV^p + VHQ^q = EVp 
QHVVp + QHQVq = Em Q . 

Eliminating \I/q, 

[?KP + E(£pp = M P (16) 

where 

^ E ^ VHQ E^aHQ QHV - (17) 

In Ref. [TtJ these expressions were also used to rewrite the orthonormality conditions for 
the solutions of Eq. (j!5|) in terms of the model functions \l/p (the solutions of Eq. ()16|) ). 
As they found however, if an appropriate choice of the P subspace is made, the usual 
orthonormality condition for ^ p can be applied. In this case the standard CI procedure can 
be used to solve Eq. ([TBI . 



D. Many-body perturbation theory 

Here we will generate a perturbation expansion for S. Define a single particle Hamiltonian 

by 

h a\\0) = e4|0) (18) 
€i = (i\ h DF \i) . 

where we introduce the operators a\ (a^) to create (annihilate) a particle. For particles in 
the core, the functions are those of Eq. (J7J). The many-body zero-order Hamiltonian is 

H = E covc + J2i a l a i} e i- ( 19 ) 

i 

where the brackets {...} denote normal ordering with respect to the closed-shell core. 
The exact Hamiltonian is 

n = J2 hnuc + E — ( 2 °) 

i i<j r *i 

where h nnc = col ■ p + ((3 — l)m e c 2 — — . H. can be separated into zero, one, and two-body 
parts: 

H m = E 

' L -^core 

= £Rt % } 



core 



(i\ h nuc \j) 



+ J2{( im \ r l2 \j m ) - ( im \ r l2 \ m j) / 
m 

(^ CI \j) (2i) 



'.i 



H {2) = 52{ala} ai a k } (ij\r£ \kl) . (22) 

ijkl 
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Expanding Eq. (|T7|) in the residual Coulomb interaction, V = H — Ho, we obtain 



E(E) = VHQ l — QHV 
t, — Al 

+PHQ—^—QVQ—^—QHV + ... (23) 

£j — rto -ft ~ "0 

One advantage of this formalism is that ho is not necessarily the same as h DF . Thus we 
may in principle use any set of functions we like in the virtual basis as long as they are 
orthogonal. In practice, however, it is important that V not get too large, and this requires 
that 

V (1) = {i\h cl - ho \j) (24) 

is small. 

We can write £ in matrix form: 



, _ „ (I\H\M)(M\H\J) 



E — E M 



+ (I\H\M)(M\V\L)(L\H\J) | (25) 

M,LeQ (E-E M )(E-E L ) 
= (E 2 ) 7J +(E 3 ) /J + ... (26) 

where I and J enumerate determinants from the subspace P, and M and L are determinants 
from the subspace Q. 

In this paper we calculate £ to second-order of the perturbation expansion. For the one- 
valence-electron case it has been shown that this level of perturbation theory, when used 
with the finite-field scaling method, is sufficient to obtain accurate results for isotope shift 



21[. Furthermore, previous studies have shown that energies calculated to this order are 



17. 1 



very accurate 

Substituting £2 into we obtain the equation of the combined CI and MBPT method, 
which we write in the same form as Eq. (|T2"j) : 

Hu + MM ) Cj = EC „ (27) 




MeQ 



Thus our method includes the core-correlation effects by simply altering the matrix elements 
in the CI calculation. In practice, this involves adding the matrix element of the sigma 
operator to the one and two-particle Coulomb integrals. 



E. Diagrammatic technique for calculating E 

Here we present all second-order diagrams for £, represented by Goldstone diagrams (see, 
e.g.: Ref. |22|). The diagrams are valid for any ion and any choice of core. Unlinked lines 
are omitted, since they correspond to states of the valence electrons that are not involved 
in the interaction. The omitted lines do affect the value of the MBPT correction via the 
Pauli principle; however, as noted in Ref. the Pauli principle can be ignored due to the 
exact cancellation of unphysical terms in different diagrams. This rule works for all orders 
of perturbation theory (this is Wick's theorem, see e.g.: Ref. |2^|). 
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The second-order diagrams can be grouped according to how many external lines (valence 
electrons) they have. There are four one- valence-electron diagrams, all shown in Fig.^ They 
correspond to the "self-energy" correction, which describes correlations between the valence 
electron and the core. Additionally, there are five so-called "subtraction diagrams" for the 
self-energy. These are diagrams that involve the external field (Eq. (|24j) ) . and are so 
named because in the Hartree-Fock field enters with a minus sign. In our formulation 
we instead use as the external field (see Eqs. ([21)1 and ([25))) which is equivalent when 
calculating £2 because VHqQ = 0. Three of the TP^ subtraction diagrams are shown in 
Fig. 12 the other two are the mirror image partners of diagramsEll and[2j2. 

FIG. 1: One-valence-electron diagrams of £ (E^ 1 )). 



a a 




FIG. 2: Subtraction diagrams ofEW. Dia grams 1 and 2 have complementary diagrams obtained 
by reflection in the vertical plane. 

X 



1 




1 2 3 



Diagrams with two external lines correspond to screening of the valence-valence inter- 
action by the core electrons. There are nine diagrams represented in Fig. |21 and four 
subtraction diagrams represented in Fig. 0J The resulting corrections to the interactions 
between valence electrons can be written in the form of an effective radial integral, as is 
usually done for the Coulomb integrals. However, it is important to note that the box dia- 
grams (Figs. E14 — EI 6) have less symmetry than the Coulomb integrals, because swapping 
the initial and final states in either the upper or lower lines changes the integral. This 
approximately doubles the number of independent radial integrals that need to be stored 
for the CI problem. In addition, for the box diagrams the multipolarity of the Coulomb 
interaction need not follow the usual rule: £(l a + l c + k)£(l b + l d + k) (see Appendix lA"j) . 
Instead, k need only satisfy the triangle conditions and can be both odd and even. This 
would again double the number of independent radial integrals, except that we have found 
that the diagrams of "wrong" parity are unimportant for carbon and may be omitted in 
order to reduce the complexity of the calculations. 



S 



FIG. 3: £( 2 ); Dia grams 1, 2 and 3 have complementary diagrams obtained by reflection in the 
vertical plane. 
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FIG. 4: £( 2 ) subtraction dia grams. Each has a complementary mirror-image partner. 
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It is worth noting also that the two-body diagrams with the largest contribution to the 
energy (the direct diagram, Fig. EJl, and its mirror) make no (linear) contribution to the 
SMS because they cancel. The exchange diagrams do contribute to the SMS; nevertheless, 
it is because of this cancellation that the contribution of E^ 2 -* to /csms is generally smaller 
than that of . 

Figure El shows a diagram with three external lines, S^ 3 \ where three valence electrons 
interact via the core. The diagrams of this type are easy to calculate (having only one 
internal summation and no summations over virtual states), but the number of corresponding 
effective radial integrals is huge. To include this diagram everywhere in the CI calculation 
would involve taking into account not only corrections to Ti.^ and Ti^, but it would also 
introduce an effective H^. Fortunately, as explained in Ref. [l7|, it is possible to omit these 
diagrams entirely provided that one makes an appropriate choice of the atomic core. We 
have not included in this paper. 
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FIG. 5: Effective three- valence-electron interaction of £2- 
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Analytical expressions for the diagrams of Figs. ^-0]are given in Appendix iBl 
A final point worth mentioning is the definition of the energy denominators in the ex- 
pressions corresponding to the diagrams. Our formalism corresponds to Brillouin-Wigner 
perturbation theory, and so there is an explicit dependence on the energy of the whole atom: 
E in Eqs. flT7j) and (J27|). 

From Eq. (|T3j) . E = E corc + E vgi , and so the core part cancels in the energy denominator 
E — Ho (see Eq. (J23j) ). The energy E ya \ corresponds to the energy of all the valence electrons. 
To calculate 7io, however, one should take into account the state of all valence electrons 
for each disconnected diagram. This is computationally impractical, but again one can 
simplify. In this paper all connected diagrams are evaluated at the energies which correspond 
to the main configuration. For example, when calculating Fig. ^1 the usual Rayleigh- 
Shrodinger perturbation theory would give the denominator (e a + e n — e a — ep); we instead 
use (e 2s + e n — e a — ep) if a and b are s-wave, or (e 2pi/2 + e n — e a — ep) if a is p-wave, and so 
on. 

The Brillouin-Wigner formalism used in our method has two major advantages over the 
Rayleigh-Shrodinger perturbation theory. Firstly, we wish to preserve the symmetry of the 
CI matrix when S is added, so we cannot have an energy denominator that depends on 
which state is initial and which is final. Secondly, Rayleigh-Shrodinger theory does not 
allow a large P space as it will produce small denominators for excited configurations. That 
is known to lead to "intruder states" - unphysical states that can lie even below the ground 
state. By contrast, Brillouin-Wigner theory is known to have the wrong limit for an infinite 
number of valence electrons. It is possible to formulate the theory without these drawbacks 
by modifying E (Ref. [23]). However, the energy dependence of the effective Hamiltonian 
is a higher order effect; compared to the prescription above, it leads to a relatively small 
energy correction that we neglect. 



IV. CALCULATIONS AND RESULTS 



We present results for each species of carbon separately, in order of increasing number of 
valence electrons. For all calculations we use a relativistic B-spline single-electron basis set 
similar to that developed by Johnson et al. ^2^, HHj 2^1 • I n all cases the MBPT calculation 



considers Is 2 as the core, and all other states as valence. We have used two different B-spline 
basis sets, each with a different number of shells included in the self-consistent Dirac-Fock 
procedure. The first set, £>i, is formed in the potential of the closed shell Is 2 core; this 
corresponds to iV DF = N core . For this set the subtraction diagrams of Figs. 121 and |U are zero. 
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The second basis set, £>2, is formed in the potential of the ls 2 2s 2 core; thus A^df > N coie 
and the subtraction diagrams must be included. 

The basis sets are described by giving the largest principal quantum number for each 
wave, e.g. 8spd6f includes the orbitals 1 — 8si/ 2 , 2 — 8pj, 3 — 8dj, and 4 — 6fj. For the 
MBPT we have used the basis AOspdfg for all ions. This basis is fully saturated in the sense 
that the energies and isotope shifts do not change with the addition of more basis functions. 



A. CIV 



ClV has one electron above a closed Is 2 core. It can therefore be treated as a single 
valence electron atom using MBPT, or as a three electron atom using CI. We have used both 
methods; the results are presented for comparison in Table In each case, the calculations 
were done using the B\ basis set. The isotope shift results have also been compared with 
previous theoretical approaches. Multiconfiguration Hartree-Fock - CI calculations were 
presented in Ref. ^3|, which also combined the Hylleraas and full-core plus correlation 
(FCPC) calculations of Refs. j2?J and [28| respectively. 



TABLE I: Comparison of energies and /csms transitions from the ground state in C IV using various 
methods. The transition energies presented do not include the addition of mass shift effects. Note 
that all results presented by other groups are non-relativistic and hence do not distinguish fine- 
structure. 

Level Energy DF + £^ Full CI Other works 
(expt.) 

Energy (cm -1 ) 
64551 64594 64564 a 
64681 64725 64399 fe 
64449 c 



2p 2 P° /2 64484 
2p 2 P° /2 64592 



&SMS (GHz-amu) 
2p 2 P° /2 -4511 -4521 -4526 a 

2p 2 P° /2 -4504 -4514 -4527 6 



-4528 



c 



a MCHF-CI: Carlsson et at, 199£ 
^MCHF: Godefroid et al, 2001 , , , 

c Hylleraas + FCPC: Results of King, 1989 j2J, and Wang et al, 1993 [2j|, combined and presented in 



Ref. 



As noted previously, the MBPT basis was completely saturated, however we have only 
included second order diagrams in our calculation. By contrast, the CI calculations are 
complete (although they do not include the Breit interaction and QED effects), but the basis 
is not completely saturated. We have used an effective 22spd8f basis for the CI calculation, 
including only single and double excitations from the leading configurations (SD-CI). We 
included triple excitations for a smaller basis, 14spd8f: this made a difference of less than 
2 cm -1 in the transition energy, and less than 1 GHz-amu for /csms- Other /-wave and higher 
partial waves were also found to be unimportant. 
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All methods give the same value for the transition energies and SMS constants to within 
0.5%. 



b. cm 



The ground state for Cm is ls 2 2s 2 We have done our calculations both as a four- 
electron CI problem (full CI) and by combining two- valence-electron CI with MBPT, consid- 
ering Is 2 as the frozen core (CI + X^ 1 ' 2 )). All results are presented in Table ITT1 Also included 
are CI results (the pure two-electron CI method) and CI + results (that do not include 
£( 2 - ) ). This allows us to examine the roles of the different parts of the core-correlation. The 
CI and CI + results are calculated with the i5 2 basis set; the complete CI + X^ 1,2 ) 
results have been calculated using both the £>i and B 2 sets. Additionally we have presented 
the MCHF results of Ref. 0. 



TABLE II: Comparison of energies, &smS) and q- values for transitions from the ground state in 
ClII using various methods. Note that the MCHF results are non-relativistic and hence do not 
distinguish fine-structure. 

Level Energy CI CI + £W CI + Full CI MCHF a 

(expt.) B 2 Bi 

Energy (cm -1 ) 
52322 52349 52383 52506 52280 
52357 52383 52418 52534 
52427 52453 52488 52592 
103365 102725 102775 103109 102598 



2s2p 3 P ° 52367 52750 
2s2p 3 Pf 52391 52784 
2s2p 3 P 2 ° 52447 52852 
2s2p l P? 102352 103719 



&SMS (GHz-amu) 

2s2p 3 P ° -3439 -3478 -3473 -3470 -3483 -3475 

2s2p 3 Pf -3438 -3476 -3472 -3468 -3480 

2s2p 3 P 2 ° -3434 -3473 -3468 -3465 -3474 

2s2p ipf -2688 -2759 -2790 -2784 -2882 -2817 

q (cm -1 ) 

2s2p 3 P$ 75 74 74 

2s2p 3 Pf 109 108 108 

2s2p 3 P$ 177 178 178 

2s2p l Pf 165 165 165 



a Jonsson et ai, 1999 Q 



For the full four-electron CI method we used a very large basis lQspdf, in the SD-CI 
approximation. This was not enough to saturate the basis entirely, and we could go no 
further because the Hamiltonian matrix became too large. The error in /csms from the full 
CI calculation could be as large as 100 GHz-amu. Nonetheless, they are in agreement with 
the calculations of Ref. Til , as well as the results of our own CI + MBPT. 



The CI + MBPT method is particularly accurate for ClII for two reasons. Firstly, 
because there are only two valence electrons there are no triple excitations, which keeps the 
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Hamiltonian small without the need for approximations (for ions with more valence electrons 
we have used the SD-CI approximation). Because the Hamiltonian is relatively small, we 
can completely saturate the basis at 20spdf. Also, there is no (corresponding to Fig. EJ); 
as stated before, in this paper we have not included it anyway. These points hold true for 
all two- valence-electron atoms; in particular we have previously shown that excellent results 
arc attainable for Mgl (Ref. [19|). 

In Ref. [01 the MCHF results were given an error of 1%; our CI + MBPT results are 
within this range, and so we believe that we have obtained a similar accuracy. It is also 
worth noting again that we have excluded the extra box diagrams with "wrong" parity from 
the results presented. The inclusion of these diagrams in S^ 2 ^ makes a difference of around 
0.1% to the /csms constants. 

C. CII 

We have treated C II as a three-valence-electron ion; its ground state is 2s 2 2p 2 P°/ 2 - We 
have used the £> 2 basis 20spdf, which corresponds to the V N ~ l potential, and we have 
restricted ourselves in the CI problem to single and double excitations from the leading 
configurations 2s 2 2p and 2s2p 2 . 

In Table IIHI we present all results for C II. Again we have presented the breakdown of 
the various parts of the CI + MBPT method. We have also performed our calculations 
using the B\ basis: this changed the results by less than 1% for all results except for the 
2 S\j2 transition, in which the difference was around 3%. For this transition, neither basis set 
was enough to completely saturate &sms- Furthermore, the difference between the results 
of CI + MBPT and MCHF-CI is fairly large for this transition (around 7%). Adding the 
next most important configuration, 2s 2 3s, to the leading set changes the energy of the 2 Si/2 
transition by 30 cm -1 (0.03%) and /csms by 14 GHz-amu (around 1%). The effect on all 
other transitions was much smaller. 

D. CI 

In Table HVl we present energies for transitions in neutral carbon. The ground state of Cl 
is 2s 2 2p 2 3 P . We used the £>2 basis of size 16spdf, and took all single and double excitations 
from several leading configurations. The energies obtained for this atom are not as good as 
those of the other ions. It is testimony to the power of the B-spline basis, however, that 
the levels appear in the correct order (with the exception of some fine structure), which is 
remarkable considering that the spectrum is very dense. 

We have generated q- values for Cl because the previous calculations (Ref. [ij|) had large 
uncertainties. We believe the new values, presented in Table II VI have errors of around 
10 cm -1 . 

In Table fVl we present the SMS constants for C I. We are within 1% of the values obtained 
using the MCHF-CI method (Refs. 0] and 0]). For most transitions the effect of core- 
correlations on /csms is around 1 or 2%, however in some cases they are larger (for example, 
in 2s2p 3 3 P° the core correlations are 8% of the total). 
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TABLE III: Comparison of energies, &sms> an d q- values for transitions from the ground state 
in ClI. Note that the MCHF-CI results are non-relativistic and hence do not distinguish fine- 
structure. 



Level 


Energy 


This work MCHF-CP CI 6 




(expt.) 


CI CI + £W CI + S^' 2 ) 



2s2p l 4 Pj 



1/2 



2s2p l % 



2s2p 2 A P, 



3/2 



2s2p 2 2 D 



5/2 



2s2p 2 2 D 



5/2 
3/2 



2s2p 2 2 S 1/2 



2s2p 2 2 Pj 



2s2p 2 2 R 



1/2 



3/2 



43003 
43025 
43054 
74930 
74933 
96494 



43118 
43144 
43186 
75587 
75585 
97095 



110624 112135 
110666 112187 



42767 

42794 

42838 

75350 

75347 

96965 

111913 

111967 



Energy (cm 
42782 
42808 
42852 
75227 
75225 
96960 
111205 
111259 



74842 

96478 
110569 



&SMS (GHz-amu) 



2s2p 2 


"l/2 


-2913 


-2956 


-2960 






2s2p 2 




-2912 


-2954 


-2958 






2s2p 2 




-2910 


-2952 


-2956 






2s2p 2 


2 D 5/2 


-2604 


-2666 


-2672 


-2672 




2s2p 2 


D 3/2 


-2604 


-2666 


-2671 






2s2p 2 


2 s 1/2 


-1204 


-1301 


-1321 


-1411 




2s2p 2 


2 Pl/2 


-1323 


-1410 


-1471 


-1531 




2s2p 2 


2p 
^3/2 


-1320 


-1407 


-1468 














q (cm -1 ) 






2s2p 2 


4 Pl/2 


132 


132 


132 






2s2p 2 


%/2 


157 


158 


158 






2s2p 2 


4p 5/2 


200 


202 


202 






2s2p 2 


2 D 5/2 


179 


181 


181 




179 (20) 


2s2p 2 


D 3/2 


176 


178 


178 




176 (20) 


2s2p 2 


2 S\/ 2 


165 


167 


168 




161 (20) 


2s2p 2 


^1/2 


162 


163 


163 






2s2p 2 


2p 
^3/2 


215 


217 


217 







"Jonsson et al., 1996 
^Berengut et al, 2004 |2 



V. CONCLUSIONS 

In this paper we have presented, in detail, a method for calculating isotope shifts and 
relativistic shifts in atomic spectra. The method uses the finite-field method to reduce the 
problem to that of an energy calculation, which is carried out using CI for the valence 
electrons combined with MBPT for the core correlations. Having previously tested the 
method in magnesium, we have now applied it to transitions in neutral carbon, ClI, ClII, 
and C IV. In all cases we have agreement with previous MCHF and MCHF-CI calculations to 
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TABLE IV: Comparison of energies and q- values for transitions from the ground state in C I using 
various methods. Note that the MCHF-CI results are non-relativistic and hence do not distinguish 
fine-structure. 



Level Energy (cm 1 ) q (cm 1 ) 

(expt.) CI CI + £« CI + S^- 2 ) MCHF-CI CI + S^' 2 ) CI 



o 2o 2 la 
2s lp Do 


O 1 /Ti A O 

21648 


oni a n 

22140 


22216 


2266b 


o 1 t cr o b 

21753 


O O 

38 




2s2p° b% 


o ' ) ~7 4 ) cr 

33735 


33529 


66264: 


o oo a n 

33240 


o o a no b 

33498 


i A r* 

146 




2s 2p6s Fq 


60333 


5980b 


60182 


( ' ( \ "1 /I /f 

60144 




-47 




2s 2p6s F" 


b(J353 


5982b 


b0202 


z^ni a 

b0164 




-24 




2s 2p6s F!} 


60393 


598bb 


b0243 


60206 




26 




o 2o o Ino 

2s 2p6s F" 


61982 


bl587 


bl975 


61911 


62002 


1 




O O 3 3 7*TiO 

2s2p° L>3 


64087 


b4773 


f* A Or\ O 

b4b28 


64562 




1/1/1 

144 


1 cr 1 / i~r\\ 

151(60) 


o o 3 3t~ho 

2s2p° °JJ" 


£? /i nnn 

64090 


b47b2 


b4bl7 


64551 


64036 


137 


1 A 1 / pn\ 

141(60) 


o „ o ~.3 3 t~\o 

2s2p D% 


a a (\(\ i 


b47bb 


A COO 


b4555 




1 /i n 

14U 


1 A cr / rn\ 

145(60) 


o o 3 3 t~>o 

2s2p F" 


75254 


76209 


7619b 


^7/3 ICO 

7bl53 




117 


111/ rn\ 

111(60) 


o o 3 3t~>o 

2s2p F% 


t cr o r cr 

75255 


76214 


76202 


7bl58 




1 O 1 

121 




o o 3 3 t~>o 

2s2p° °Fq 


7525b 


76207 


76194 


7bl51 




1 1 cr 

115 




2s 2p6d L>2 


77680 


79297 


79724 


79b43 




i 




Is 2p4:S Jtq 


1 n CT 

78105 


79737 


80152 


8007b 




no 

-33 




o 2 o /( 3no 

2s 2p4s ±i 


78117 


79750 


oni r* o 

8016b 


o c\c\c\r\ 

80090 




O 1 

-21 




o 2 o /( 3no 

2s 2pis ±2 


78148 


79787 


80205 


om on 

80130 




24 




o 2 o o j 3 Trio 

2s 2p6d J J^2 


f701 C\C\ 

78199 


^n o /i cr 

79845 


80271 


oni no 

80193 




O 1 

-31 




ZS zpou x*3 


7891 P. 
1 oZIO 




ouzoy 


8091 1 
oUZl 1 




i a 

-10 




2s 2 2p3d 3 Df 


78293 


79937 


80354 


80275 




-13 




2s 2 2p3d 3 D% 


78308 


79954 


80371 


80293 




13 




2s 2 2p3d 3 D% 


78318 


79966 


80385 


8030b 




29 




2s 2 2pAs 1 J Pf 


78340 


79983 


80403 


80327 




17 




2s 2 2p3d ^ 


78530 


80802 


80626 


80549 




15 




2s 2 2p3d ^Pf 


78731 


80402 


80822 


80746 




12 




2s 2 2p3d 3 P 2 ° 


79311 


80957 


81307 


81225 




18 




2s 2 2p3d 3 Pf 


79319 


80967 


81318 


81237 




30 




2s 2 2p3d 3 P ° 


79323 


80972 


81323 


81242 




35 




2s 2 2pM x Bl 


83498 


86521 


86942 


86858 




8 





a Berengut et a/., 2004 Vm 
^Carlsson et at, 
c J6nsson et al, 1996 [la| 

within a few percent. In Table IVll we compare our calculations with the few experiments that 
exist for carbon ions; in all cases agreement is within around 0.005 cm -1 , which corresponds 
to an error in /csms of around 20 GHz-amu. 

In Table IVlTl we present total isotope shifts for some important transitions. These tran- 
sitions can be observed in quasar absorption spectra, and can therefore be used to probe 
variation of a and isotope abundance evolution. The results are presented both in MHz and 
km/sec: the latter is the preferred form for use in astronomy. 
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TABLE V: Comparison of £;sms f° r transitions from the ground state in C I using various methods. 
Note that the MCHF-CI results are non-relativistic and hence do not distinguish fine-structure. 



Level 


Energy 




&SMS (GHz-amu) 






(expt.) 




CI + EW CI + S^ 1 - 2 ) MCHF-CI 


Zs zp Do 


Zlb4o 


186 


180 


191 


152 a 


zszp D 2 


qq7q r 
OO/ OO 


ok An 


-2579 


-2588 


-2583 a 


ZS ZpoS ±q 


bUooo 


1 ACiK 
14Ub 


1405 


1419 




ZS ZpoS ij 


bUooo 


14Ub 


1406 


1420 




zs zpos ij 


ouoyo 


1 /I HQ 

14Uo 


1408 


1422 




ZS ZpdS Tj 


blyoz 


1549 


1551 


1559 


1553 a 


n n n^3 3no 
ZsZp JJ3 


b408/ 


-zlb5 


-2224 


-2227 




n n n^3 3no 
ZsZp L>i 


b4(jyu 


ni /;c^ 

-zlo5 


-2224 


-2227 


-2222 & 


ZsZp L>2 


b4Uyl 


-zlo5 


-2224 


-2227 




ZsZp 


f7rnr /I 

/5z54 


-Iz < z 


-1390 


-1392 




o^o^3 3r>o 
ZSZp J"2 


(bZOO 


1 070 


-1389 


-1391 




o^o^3 3no 
zszp Jq 


tuZoo 


1 971 

-IzYl 


-1390 


-1392 




zs zpoa 


1 fboU 


loo4 


1320 


1331 




ZS zp4s Xq 


70 1 n^ 


1 QOQ 

lays 


1392 


1407 




ZS Zp'i.S ij 


/Oil / 


14U4 


1397 


1412 




9 ^2o/v^/( ^ 3 do 
ZS Z£>4S fn 


7Q 1 y|Q 

(8148 


i /) i £ 

141b 


1408 


1422 




Zs Zp6a i<2 




1 Q7Q 

13 (o 


1368 


1381 




Zs Zp6a r 3 


7001 £> 

/8zlb 


1381 


1372 


1384 




o^2o/vaQ^7 3no 
zs zp3a i/J 


/ozy3 


1430 


1 A 99 
14ZZ 


1 A Q/1 
14o4 




o^2o/vaQ^7 3no 

zs zp3a 


/8308 


i A on 
14zy 


1421 


1432 




2s 2 2p3d 3 Df 


78318 


1430 


1420 


1432 




2s 2 2p4s 


78340 


1443 


1435 


1446 




2s 2 2p3d ^ 


78530 


1451 


1440 


1452 




2s 2 2p3d 


78731 


1436 


1426 


1438 




2s 2 2p3d 3 P 2 ° 


79311 


948 


998 


1010 




2s 2 2p3d 3 Pf 


79319 


956 


1006 


1018 




2s 2 2p3(i 3 P ° 


79323 


960 


1009 


1021 




2s 2 2p4ci ^ 


83498 


1277 


1258 


1268 





"Carlsson et al. 
kjonsson et al., 



1995 
1995 
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TABLE VI: Comparison of calculated C - C isotope shifts with experiment. 

Transition A 5^ 13,12 (cm -1 ) 

Lower Level Upper Level (A) Expt. This work 
~CI 

2s 2 2p 2 3 P 2 2s2p 3 5 5f 2967 0.670(5) a 0.674 6 
2s 2 2p 2 % 2s 2 2p3s tyf 2479 -0. 156(3) c -0.151 

-0.156(2) d 

CII 

2s2p 2 2 S 1/2 2s 2 3p 2 P.° /2 2837 -0.612(2) c -0.617 
2s2p 2 2 S 1/2 2s 2 3p 2 P° /2 2838 -0.623(3) c -0.617 

"Bernheim and Kittrell, 1980 

b Actually, the 2sl2p 2 3 P Q - 2s2p 3 5 S% transition was calculated. 
c Burnett, 1950 £3 
d Holmes, 1951 jij 



TABLE VII: Total 13 C - 12 C and 14 C - 12 C isotope shifts for important transitions. We believe 
the errors are of the order of 0.1 GHz. 

Transition A 5^ 13 ' 12 5u 14 ' 12 

Ground State Upper Level (A) (GHz) (km/sec) a (GHz) (km/sec) a 



2s 2 2p 2 3 P 


2s 2 2p3s 3 Pf 


1657 


-2.75 


0.46 


-5.09 


0.84 




2s2p 3 3 Df 


1560 


21.10 


-3.29 


39.12 


-6.10 




2s2p 3 3 Pf 


1329 


16.91 


-2.25 


31.34 


-4.17 




2s 2 2p4s 3 Pf 


1280 


-0.82 


0.10 


-1.51 


0.19 




2s 2 2p3d 3 D^ 


1277 


-0.94 


0.12 


-1.75 


0.22 




2s 2 2p4s l P? 


1276 


-1.01 


0.13 


-1.88 


0.24 




2s 2 2p3d 3 Pf 


1261 


1.84 


-0.23 


3.42 


-0.43 


CII 














2s 2 2p 2 P° /2 


2s2p 2 2 D 3/2 


1336 


25.10 


-3.35 


46.53 


-6.21 


2s2p 2 2 D 6/2 


1336 


25.10 


-3.35 


46.54 


-6.21 




2s2p 2 2 S l/2 


1037 


18.70 


-1.94 


34.66 


-3.59 


cm 














2s 2 % 


2s2p l P{ 


977 


28.76 


-2.81 


53.33 


-5.21 


CIV 














2s 2 S 1/2 


2p 2 P° /2 


1551 


35.89 


-5.57 


66.54 


-10.32 




2p 2 P° /2 


1548 


35.79 


-5.54 


66.35 


-10.27 



a 5v = 5X/X x c (km/sec). A negative velocity means that C (and C) are at lower wavelength than C. 
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APPENDIX A: MATRIX ELEMENT OF THE TWO-BODY OPERATOR 

The two-body operator used is this work is the sum of the Coulomb interaction operator 
and the "rescaled" SMS operator (atomic units): 

Q = ] r + Api p 2 = J2Qk, ( Al ) 

|ri-r 2 | k 

where A is the scaling factor, p = — zV is electron momentum, and 

Qk = — T -^y t (ni)n(n 2 ) + A pi • p 2 S kl . (A2) 
We use the following form for the single-electron wave function 



Here a = 1/137.036 is the fine structure constant, and fl(n)ji m = —(a ■ n)Q(n)ji m . 
The matrix element of operator (jA2|) with wave functions (jA3|) has the form 




(A3) 



(0a(ri)0 6 (r 2 )| Qk \M r l)M T 2)) = Cab,cd(Rab,cd ~ ^klPacPbd), 



(A4) 



where the angular factor Ck is the same for both operators 




(A5) 




1, if x is even 
0, if x is odd 



'ab,cd 



is radial Coulomb integral 



ab,cd 




(fa(ri)f c (n) + a 2 9«(ri)gc(ri)) (Mr 2 )fd(r 2 ) + a 2 9b(r2)9d(r 2 )) dr x dr- 



2, 



(A6) 



while pi2 is the radial matrix element of the SMS operator 



P12 



(A7) 
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APPENDIX B: EFFECTIVE INTEGRALS FOR S DIAGRAMS 



Here we give expressions for the diagrams in Figs. ^ 0] The following conventions for 
the indices are used: a, b, c, and d correspond to the external lines (valence electrons); 
1 < m, n < iV corc ; N corc + 1 < a, (3. Following 22|| we use the shorthand notation [j, k, 



(2j + l)(2k + 1) . . .. Note also that all Wigner three-j symbols come with a parity selection 
rule, which is not explicitly given. This is because the three-j terms come from the formula 

Additionally, the subtraction diagrams (Figs. 121 and 0J) have terms which are matrix elements 
of h CI (see Eq. ()14|) ). This operator cannot change the angular momentum of the state; for 
example (a\ h cl \n) comes with a factor 5j a j n 8i a i n which is not presented explicitly. 

We start with the contributions to the self-energy matrix elements from the diagrams of 
Figs. ^ and EJ All of these expressions have the following additional terms that we do not 
write explicitly: 5 jajb , 5 mamb , and 5 lah . 

Fig. [Hi: 

n \ - [jnjajp] ( j a jf3 k\ 2 (j n j a k \ ^ Rgn^a R Pa,bn / R ,x 

D «* = <L Fli i/o _i/o n i/o _i/o n — TZ ; ~ ^l) 



nap 
fc 



Fig.[U2: 

Ua > b ~ ^ ' v»>3°>m { l/ 2 -l/ 2 M 1/2 -1/2 j I 1/2 -1/2 



2 



na/3 



v f ja ja k 2 \ f ja j/3 h \ Rgn,p Q R 0a,nb / R9 n 

I X h -V* ) \ j n j a hje a + e n -e a - e p { ^ ] 



Fig. 033: 



p. [jmi jni ja] ( ja jm k\ f jn ja k\ R aa,mn R mn,ba fRQ\ 

a ' b ~ ^ [k] [ V2 -y 2 O) [i/2 -i/2 e m + e n -e Q -e b {B6) 



run a 
k 



Fig. [U4: 

7-j _ t_-\ Nfci+fca+lr„" „• „• 1 | j™ Kl| [ Ja Jn | / Jm Ja &5 

' 6 ~ ^ [Jm,Jn,Jaj I y 2 _y 2 q I I l/ 2 _l/ 2 Q J I l/ 2 _l/ 2 Q 



fc 1 fc 2 



jn Ja fe 2 \ / Ja jn &2 1 R aa,mn R 



fe 2 

mn,«6 



Fig. 01: 



Fig. |2l2: 



w Jn Ja 'T2 I Ja Jn '"2 I "tta,inn--mii,aii (B41 

^ 1/2 -1/2 j \ j Q j m fci J e m + e n - e Q - e b 

A*,6 = 2^«] '— — (B5) 



D« = f iy 2 V (B6) 



fc 
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Fig. H3: 



^ = jWW^) (B7) 



Next we present contributions to the valence-valence screening diagrams due to the dia- 
grams of Figs. and 0J These diagrams have fc as an external parameter (not to be summed 
over). In the CI procedure they are added to R!l bcd of Eq. (|A4|) . 

Fig. 01: 

jjk b«'^™] ( jot jn k\ R a n,ca^ab,nd /X3Q\ 

Kb,cd - 2^ ri.i i/o -i/o n — ; — ; — ^ 



[fc] V Va -1/2 

Fig. 02: 



= £(-i) fcl+ U,-7»] 



ja jn fcl W ja jc Mf jn ja fc 
1/2 -1/2 ill/2 -1/2 ill/2 -1/2 



x( h H, n ) I 3a 3c 1 Rk ^ ncRknbM (B9) 



Fig. 03: 



R k ab , cd = E(- 1 ) fcl+fe b-^] 

na 



jb jn fcl W ja fcl W jn ja fc 

1/2 -1/2 ill/2 -1/2 ill/2 -1/2 



v ( 3b 3d k\ 1 j J b j d k \ Rto,nd R na,ac fvnn \ 
HV'-V'Oj (j a j n fcl | e „- ea (B10) 

Fig. 04: 

-ic+id+ja+jn [„■ „' LI f ja jn fcl | | j" j<i fcl 1 
Ua,3n,K\ I y 2 _y 2 q I I y 2 _y 2 g \ 



^abfid 


E(-!)^ 

not 
k l k 2 






Jn J. 

v V 2 - ] 




H 


f jc 3a 
[ fcl fc2 


05: 






r>k 

^abfCd 


E(-i) Ja 

na 




x| 


( ja j. 

v V 2 - ] 




x < 


1 3c 3a 
[ fcl fc2 



c k 2 \ ( j b j a k 2 \ ( j a j c k \ 1 / j b jd k \ 1 

l /2 J \ 1/2 -1/2 ) \ 1/2 -1/2 y 1 ^1/2 -1/2 y 1 

fc ) f jb jd fc 1 ^ao.nd^X" 
Jn M fc l fc2 J« J e « - 



+jb+jc+jd+ic+inr„- „■ jLl f Ja Ja fcl I f J" Jd fcl 1 

L7a,J„,«] ^y 2 _y 2 g J ^ y 2 _y 2 g J 

c fc2 \ / jb jn k 2 \ ( j a jc ^ V Jfo Jd fc 

L/ 2 ^ 1/2 -1/2 J \ 1/2 -1/2 j ^ 1/2 -1/2 

fc I f jb Jd fc 1 R ln,ad R ab,cn 

3a J I fci fc 2 Jn J e n - e Q 



(Bll) 



(B12) 
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Fig. H6: 

T)k _ \^ ( _T\ja+jb+jc+jd+jm+jn+k+k 1 +k2+l\- • LA I 3 a jm fcl \ I 3 b 3n 
n ab,cd ~ L ) [3m,3n,K\ I y 2 _y 2 q I _y 2 Q 



mn 
k 1 k 2 



X 



Fig. 



Fig. |U2: 



jm 3c h \ ( j n U k 2 \ ( ja jc M f j& jd 

v 1/2 -1/2 y ^1/2 -1/2 J ^ 1/2 -1/2 J ^1/2 -1/2 y 

x f Jc Ja ) f Jd Jfe 1 Rgb,mnRmn,cd (B13) 

\ h k 2 j m j 1 h k 2 3n I e m + e n -e c -e d 



■ „ R k ancd {b\h cl \n) 

R k ab,c d = ~ E ' (B15) 
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